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A new scattering approach for correlated one-dimensional 
systems is developed. The adiabatic contact to charge reser- 
voirs is encoded in time-dependent boundary conditions. The 
conductance matrix for an arbitrary gated wire, respecting 
charge conservation, is expressed through a dynamic scatter- 
ing matrix. It is shown that the dc conductance is equal to 
e 2 jh for any model with conserved total left- and right-moving 
charges. The ac conductance matrix is explicitly computated 
for the interacting Tomonaga-Luttinger model (TLL). 

Pioneered by Landauerjil the scattering approach for 
quantum transport has proven powerful in mesoscopic 
physics. Nevertheless, it is restricted to non-interacting 
systems, and to the stationary regime. There were for- 
mal extensions to finite frequency transport based on a 
self-consistent approach,!] or non-equilibrium techniques 
for interacting dotsja but these formalisms are difficult 
to exploit. Proposed here is a new scattering approach 
at arbitrary frequency for linear transport through a 
strongly correlated, one-dimensional wire in the low- 
energy regime. Charge reservoirs connected adiabatically 
to the wire are accounted for by appropriate boundary 
conditions. Coupling to a gate is taken into account, 
ensuring charge conservation. The corresponding AC 
3x3 conductance matrix is expressed through a novel 
dynamic "scattering" matrix S(w). Further progress is 
then made in two cases. First, for any model where 
the total charge for right- and left-moving electrons is 
conserved, the transmission is shown to be unity in the 
zero-frequency limit. This generalizes the DC conduc- 
tance result a. =j= e 2 jh shown for a Tomonaga-Luttingex 
liquid (TLL)ol3 or arbitrary finite-range interactions.Ll 
The same result was shown in Ref. [ [| through different 
arguments restricted to the stationary regime, without 
describing the reservoir- wire interface. Second, S(u>) is 
computed for the TLL model, giving an AC conductance 
that depends on interactions in contrast to the stationary 
regime. 

Without connecting one-dimensional leads to an inter- 
acting wire, this work extends the concept introduced 
m Refs. [ where reservoirs are simulated by the 

electrons they inject. The leads have served to define 
the incident and transmitted electrons, different from the 
proper modes of the wire. S(w) will be related to the dy- 
namic transmission. For a TLL model, the same conduc- 
tance results were found by the author by computing the 
curren±-|in response to an appropriate external electric 
field, PU 1C I More recently, they were confirmed by Blanter 
et aLJ through a self-consistent treatment of interactions, 
justified in the absence of backscattering. Other works 
based on the Kubo formula in a TLL with leads found 



different results due to a different electric field profile to 
which current is very sensitive. tn 

An underlying hypothesis of Landauer's approach for 
noninteracting systemscl is the ideal nature of the con- 
tacts, ensuring that emerging electrons are absorbed 
without reflection by the reservoirs .u Such a concept can- 
not be extended to interacting systems, as illustrated 
in Ref. [ Q.Ejo Rather, interactions give rise to collec- 
tive excitations, or Laughlin quasiparticles in edge states, 
that are different from the electrons in the reservoirs. 
An emerging 'jouasiparticle" undergoes a quasi- Andreev 
type reflectioroCi'EiJ'lij at a perfect contact with a reser- 
voir. 

This paper is mainly concerned with systems con- 
nected locally to reservoirs, such as quantum wires, 
nanotubes; edge states couple differently to reservoirs 
and the present scattering approach will be extended 
elsewhere. EJ 

Consider an arbitrary one-dimensional finite wire de- 
limited by [—a, a]. The long wavelength part of the elec- 
tronic density can be decomposed into right and left- 
moving electron densitiesllll p + and p_ including implic- 
itly the zero modes, p = p + + p_ where spin is ig- 
nored for simplicity. For r = ±, the boson field $ r 
defined by p r = — d x & r /2Tt is the canonical conjugate 
to rp r (Kac-Moody algebra). The kinetic Hamiltonian 
is Hkin — J_ hvp(p+ + /0_)/2. Any interaction Hamil- 
tonian Hi n t either between electrons or with impurities 
can be expressed as a functional of <£>_, thus the total 
Hamiltonian. 



H = H kln + H int + eVgateQ = H(p + , p-), 



(1) 



is a functional of p±. Coupling to a gate is incorporated, 
and Q = Jl a p(x). 

The current field j(x) can be expressed independently 
of the dynamics, in or out-of-equilibrium. For this, the 
Hamiltonian H^ in the presence of a vector potential 
A, is used : 



dA(x) 



(2) 



A=0 



A can be absorbed by a gauge transformation of the right 
and left-going Fermion fields, 'Jv ~ e lr t[ for r = ±.EZl 
This is accomplished by the substitutionEl 

$ r (a?) -> * r (:c) - — / A{x')dx'. 

Taking the spatial derivative, one obtains H^ as 
H^=H( P++ e -A,p_- e -A 
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Differentiating with respect to A [Eq. (§)] yieldsEl 
For r = ±, p r are operators that play a central role 



(3) 



Also of use will be their sum, 



V 3 (co) 



(5) 





FIG. 1. A gated wire connected adiabatically to charge 
reservoirs with time-dependent electrochemical potential 
£ti,2(i) = e 8 "*^ii,2(w). The boundary conditions apply in the 
presence of arbitrary backscattering inside the wire. 



Consider now a typical transport measurement, where 
one connects the wire adiabatically at ±a to charge reser- 
voirs [see Figure] . Spatial and temporal structure on the 
scale of \f (Friedel oscillations) are ignored. The lo- 
cal Hamiltonian density 7i(±a) is quadratic in p± when 
backscattering is irrelevant, so that p± (±a) are linear 
in p± . One can think of the expectation value of p+ (x) 
as the energy required to add a right-going electron at x. 
It includes the local Fermi energy hvFP+(x), the interac- 
tion energy, and the gate potential that shifts the bottom 
band. In some sense, it is a local electro-chemical poten- 
tial for right-going electrons. On the other hand, the left 
(right) reservoir injects bare right (left)-going electrons 
with a well defined electro-chemical potential p\ {p2)- 
At the contacts where both incident fluxes impinge, the 
energy is conserved in the absence of any dissipation pro- 
cesses. Thus the field p + (—a,t) [p^(a,t)} is required to 
be pinned to p\ [^2] at any time. To extend this con- 
jecture to alternative regimes, time variation has to be 
slow enough so that the reservoirs are driven adiabati- 
cally through a sequence of equilibrium states with well 
defined time-dependent electro-chemical potential: 



(i + (-a,t) = m(t) 
P-(a,t) = ii 2 (t). 



(6) 



Without interactions near the contacts, p±{ z ^a,t) — 
hvFP±(T a ,t) + eV ga te{t)\ Eq. (||) imposes the density 
for incident electrons, and generalizes the Landauer con- 
cept to AC transport with arbitrary backscattering. But 
in the presence of interactions, p± depend on both p + 
and p-, thus the density for incident electrons is not im- 
posed, in contradiction with Ref.[ Rather, electrons 



are partially reflected, giving rise to a contact resistance. 
Indeed, Eq. (H) leads to a discontinuous local electro- 
static potential Vi oc . eVi oc follows the electro-chemical 
potential on the reservoir side, while 

\r t \ 9[H-H kin ) hvp 
eVi oc (x) = Q—r = fi(x) - —P{x) (7) 



on the interacting side.EJuNote that eV ga te is included. 
p(x) is. a local electro-chemical potential for both 
carriers!] [Eq. (0)]. This clarifies a confusing point in 
the self-consistent treatment of interactions in higher- 
dimensional systems, where the continuity of Vi oc is ex- 
pected, even though not implemented in the results. It 
would be interesting to define and then verify analogous 
conditions to Eq. (Q). 

The AC conductance matrix Gs(w) will be now ex- 
pressed formally, where u> is the external frequency. In 
the sequel, I switch to the Fourier transform of expecta- 
tion values in a time-dependent ground state. Gs(w) is 
a 3 x 3 matrix with g a p = 5I a /5V(3 where Vi,2 = Pia/ e i 
V3 = Vgate, Ii,2 = TjiTa-i^) an d ^3 the displacement 
current. A complete description with all the surrounding 
three-dimensional environment would be too complex. 
Instead, assume that all the electric field lines emerging 
from the wire end up on the gate, thus the latter carries 
an opposite charge to that on the wire Q(ui) . This ensures 
Kirchoff 's law, I a — because, using the continuity 
equation, 



Iz{u>) = —ieLoQ{u}) — I iujep{x,uj) (8) 

J —a 

d x j(x,u>) = -h(u>) - 



V ga te{u) appears as a reference po- 
, Thus the two constraints on the 



On the other hand, 
tential in Eqs. 
conductance matri: 



9*0 = = ^ 9af: 



(9) 



are ensured. Next focus on the first 2x2 block of Gs{co) 
denoted G 2 (w). Using Eqs. (§,|), 

h = -j(-a,Lu) = - [/x_(-a,w) - MiM] , 

h = j{+a,uj) = - [p + (-a,uj) - p, 2 (u)} . (10) 

In order to express I a to linear order in p,± t 2 = 
p,±(^fa,uj), it is sufficient to retain the linear dependence 
of /z±(±a,Ci;), determined by some matrix S(c<j), 



H+{a,uj) J \ p-{a,u) 1 



Combined with Eq. (Oil), this gives 



g 2 m = -[sm-i] 



(12) 
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In order to interpret S(u>) as a "scattering" matrix, 
one can model the transition region by an N-chamiel 
system where electrons are injected at the same quasi- 
cquilibrium distribution. On the reservoir side of —a, 
the total current is J+(w) — J_(o;), where J+(lu) = 
Nefii(u)/h is the incident current, and J—(u>) the un- 
known reflected current. The local continuity of the 
current at —a together with Eq. imply that 

efj,-(—a,oj)/h = J-(ui)/N is the average reflected cur- 
rent per channel. Similarly, epL + (a,u))/h is the average 
transmitted current to the right reservoir. If the elements 
of S(ui) are denoted as follows 



R(uj) T'(cu) 
T(ui) R'(u) 



(13) 



then T(u>) [R(u)] can be viewed as the total dynamic 
transmission [reflection] coefficient for the incident flux 
from the left to the right reservoir [into the left reservoir] . 
T" and R' play the same role for the right reservoir. Nev- 
ertheless, important differences from the usual scattering 
approach should be stressed. The elements of S(w) deter- 
mine directly the current or density, but are nonetheless 
complex numbers. S(w) is not unitary, and in general not 
symmetric unless there is a perfect reflection symmetry. 
In addition, T(lu) + R(w) ^ 1, and current conservation 
is ensured by the gate [Eq. (g)]. Finally, using Eq. (12), 
and letting the lo dependence be implicit : 



R - 1 V l-R-T' 

G :s = — | T R' - 1 l-R'-T 

l-R-T l-R'-T' R + R' + T + T'-2 

(14) 

At zero frequency, 5(0) becomes real symmetric, and 
T(0) + R(0) = T'(0) + R'(0) = 1, but R(0) can be nega- 
tive. Then the DC conductance is 



5i2 = -.9ii = T(0)- 



(15) 



As a first application of these boundary conditions, 
consider now a model where both total charges Q± = 
f p±(x)dx are conserved, 



[Q±,H] = 0. 



(16) 



Then it is shown here that T(0) = 1. For a quadratic 
Hamiltonian, this result was shown in Refs.[ §|7|, and 
justifies the hypothesis adopted in Refs.[ ||. 

In the Heisenberg representation, an operator O 
evolves according to ihdO/dt = [H, O] + ihdO/dt, but 
dO/dt = in the stationary regime. For r — ±, <I> r is 
the canonical conjugate to p ri thus d& r /dt = —rp r /h [ 
Eq. (§]. Then hdp r /dt = rd x \i r is an equation for field 
operators that one can integrate between —a and a to 
get, using Eq. (M), 



p r (a,t) - p. 



(17) 



On the other hand, the field a, t) cannot fluctu- 

ate but is equal to pi [Eq. @], thus p, + {a,t) = pi at 
all times. Similarly, a,t) — U-(a, t) = \i2- Thus 

T(0) = 1 and R(0) = [sec Eq. Q)] and the DC con- 
ductance is equal to e 2 /h [see Eq. ( ^5| ) or simply Eq. 

(!)]■ 

A second application is to investigate ac transport in 
the simplest model (verifying Eq. (|lq)) : the TLL model. 
The matrix S(u>) in Eq. ( |ll| ) can be computed in an 
instructive way, through a "transfer" matrix A(a>) such 
that 



fi(a,u>) — A(ll>)/jl(— a, u), 



(18) 



where fi stands for the vector For this, it is 

convenient to use the right- and left-propagating cur- 
rent modes j±, corresponding to up and down edge- 
excitations in a Hall bar,El that can be denoted "quasi- 
particles". The Hamiltonian density is given by 7i(x) = 
h (j^. + j 2 _^jJ2e 2 uK , where u and K are interaction 
parametersilll Without interactions, j± — evFP±, u = vp 
and K = 1. j± propagate freely at the sound velocity u, 
thus 



j(o,w) 



_ i<J z ujt L 



j(-o,<*>), 



(19) 



where = 2a /u is the transit time of the wire and <x z 
the z Pauli matrix. On the other hand, j± are related to 
p± by simple diagonalization, but their relation to p± is 
of more use here : 



M 



'- M/x(x, u) 



1 



1 - 7 
- 7 1 



1 + 7 

where the coefficient 7 is given b 



7 = 



1 - K 
1 + K' 



Then M 1 can be obtained from M by 7 
tions (|20|,[T9j) yield the "transfer" matrix 



AH 



M~ 1 e ia '*" tl 'M. 



(20) 

(21) 
-7. Equa- 

(22) 



This allows to deduce the scattering matrix S(w) in Eq. 
(|ll|), symmetric due to the reflection symmetry, 



T( W ) =!» = (! - 7 ). 



7e" 



e" 

R{lo) = R\lo) = 7 [1 - e lutL T(iu)] 



(23) 
(24) 



One can show that |DefS(ci;)| = 1, which is a constraint 
that has to hold for any quadratic Hamiltonian with time- 
reversal symmetry. Note that S(uj) depends solely on the 
intrinsic properties of the TLL model; the boundary con- 
ditions (0) allow to express the AC conductance matrix 
through S(w), Eq. (pT[). I now analyze in more details 
the capacitive effects. The gate conductance is 
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ff33 (a;) = 2-(l- 7 ) IT -^ 



(25) 



Thus the "electro-chemical" capacitance of the wire per 
unit length with respect to the gatdl2lO 



C = - lim 



533 M 



iui2a 



Ke 2 
u h 



is proportional to the compressibility.0 This result can 
be checked by minimizing the zero mode contribution to 
H at fixed total current, thus by minimizing huQ 2 /8aK+ 
eQVgate- C results from two capacitors in series: its value 
without interactions Cq — e 2 dn/dE = e 2 hvp/2, of purely 
kinetic origin, and the "electrostatic" capacitance c, ob- 
tained by evaluating Eq. (0),0Oy 

Sp 



e 2 


(-- 


h 


\K 


1 


1 


- -f 

c 


"Co' 



(26) 
(27) 



SVloc 

1 
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Interestingly, evaluating then differentiating Eqs. (|7[|5|) 
with respect to p allows to recover Eq. (|2^) . In the TLL 
p{x) = e 2 p{x)/C + eVgate justifying its interpretation as 
a local electrochemical potential for both carriers in Ref . [ 
M. But measuring C gives the ratio K/u, leaving both 
u and K unknown. For usual ballistic quantum wires, 
where a is several pm, 2it/ti, ~ GHz is quite high. One 
is oftenjn the regime cot l *C 1 , where the non-dissipative 
parMtll of 333 (w) [Eq. ©] is -Im[g 33 (uj)} ~ 2X - 
X 3 (l - l/3K 2 )/2, with X = Cato. A strategy consists 
in measuring the leading term 2X, then the subleading 
term that one divides by X 3 /2 to infer 1 — 1/3K 2 , thus 
K . u can be then determined from C . 

The underlying dynamics are now interpreted. Eq. 
(p0|) is equivalent to 



ep—(x, uo)/h 



7 1 + 7 
1-7 -7 



ep+(x,uj)/h 



(28) 



so that the matrix on the right hand side can be viewed 



ff 



as a local "scattering" matrix.Hu Let us focus for in- 
stance on x = —a where p+(— a,u>) — pi(to). When 
no charge is incident from the left reservoir, i.e. p\ = 0, 
then j + = — 7 j_; —7 is the reflection coefficient for a 
"quasiparticle" incident on the contact. For repulsive 
interactions, K < 1, thus —7 < 0; a "quasinhpic^, is re- 
flected, in analogy with Andreev reflectionO'Bl 10 H 15 l If no 
"quasiparticle" comes from the right, i.e. j_ = 0, then 
one finds 



3+ = t(1 - 7)Mi> 



(29) 



and thus K a = 2K/(\ + K) — 1 — 7 is the transmission 
coefficient for the incident flux from the reservoir. T(lo) 
and R{u>) result from the multiple reflections on-ifcbe con- 
tacts, in analogy with a Fabry-Perot resonatorBE3 They 



have resonances at the collective modes of the finite wire 
Lj n = uq for q = 2mr/2a, at which T = 1 and R = 0. 
This is because e iuJn(T ^i = I, thus Eq. (|2|) becomes 
A(u! n ) = MM -1 = l£3 Note that since j = j+ - 
Eq. (129) obtained for j_ = yields the current at the 
interface of a semi-infinite TLL and a F-Mpi liquid, the 
DC conductance becomes g a = K a e 2 /h.aa^-i 

Indeed, all the above scattering matrices have been 
encountered in Refs.[ p|,ppC|], where a TLL is connected 
perfectly to noninteracting leads at ±a. 

More generally, consider the case where Hi nt is an ar- 
bitrary quadratic functional of p, vanishing for \x\ > a. 
If an electron impinges at t — on —a, i.e. p + (x,t — 
0) = §{x + a), the transmitted [respectively reflected] 
charge to a [at —a] at time t, i.e. p+(a,t) [p~(— a,t)] 

is given by the function M ++ (a, —a, t) [M |_(— a, —a, f)] 

whose Fourier transform coincide exactly with T(u>) (re- 
spectively R(u>)) in Eq. (p8h. In addition, these func- 
tions deterirjijie the non local dynamic conductivity at 
the contacts,™] 



a(a, -a,oj) = -tT{uj) 
h 

e 2 

a(a,a,uj) = — [1 - R(u))]. 



(30) 



The reservoirs can be modeled kv— an external potential 
that drops only at the contacts,BtI£l so that, taking into 
account the constant gate potential in [—a, a], 

E(x, lo) = [Vi - V ga te} 6(x + a) - [V 2 - Vg a te} 6(x - a). 

Then j(±a,u) = a(a, =pa, u) [14,2(w) - V ga te (w]L lea 
to the same Gs(w) [Eq. (p4|)] by use of Eqs. (J30|, 

This is not a pure coincidence. The action being 
quadratic, the ground state properties are given exactly 
by minimizing it. ThepCquation of motion thus obtained 
imposes the continuityul3 of both j = e [p + — pJ\ /h and 
that of Eq. (|^), equals for instance p(x) — hup/K + 
eVgate in the TLL. Here the to dependence is implicit. 
Thus both p + and /i_ are continuous. For any |a;| > a, 
p±(x) = hvFP±(x) (see Eq. (^)). The left reservoir in- 
jects electrons with density p+{— aS~\uS) — pi(u)/hvF 
on the noninteracting lead side, thus p+(—a^\uj) = 
Pi(uj). This fixes the continuous field /z+ on the inter- 
acting side, p+(— , uj) = pi, which is exactly Eq. <^). 
Similar reasoning holds symmetrically for the right reser- 
voir. Also p-/h has to be continuous at —a, thus it is 
equal to vpp^, the reflected current. This is in accor- 
dance with the previous general interpretation of S(vX 
since one-dimensional leads model an effective channel.E-a 

These arguments in favor of the noninteracting leads 
are restricted to a quadratic density Hamiltonian H.(x) 
along the wire, although some of them might be extended 
if only H(±a) is quadratic. JLhpj-ole of backscattering in 
a TLL connected to leadsaHllJ was found to be con- 
trolled by the Fabry-Perot dynamics recovered here by 
using Eq. (^|), but the equivalence has to be checked 
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and might be limited to linear transport. The bound- 
ary conditions (0), formulated without connecting leads, 
are more general, and offer possibilities for future stud- 
ies. One has to reformulate the bosonisation procedure 
to compute the correlation functions. Implementing Eq. 
@ in a path integral formalism would give access to the 
nonlinear regime, AC transport, and current fluctuations. 
Conceptually, the scattering approach presented here can 
be extended to situations where linear response theory 
fails, as will be the subject of aJirture study on edge 
states and many channel systems. t£l 
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Note : In a recent erratumj^j Egger and Grabert mod- 
ified their boundary conditions for the TLLEl by using 
self-consistent arguments^ Their corrected results agree 
with Ref. [ Q and therefore with its present generaliza- 
tion. 
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